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ABSTRACT 



We report measurements of a polar coronal hole during the recent solar min- 
imum using the Extreme Ultraviolet Imaging Spectrometer on Hinode. Five ob- 
servations are analyzed that span the polar coronal hole from the central meridian 
to the boundary with the quiet Sun corona. We study the observations above 
the solar limb in the height range of 1.03 — 1.20 R & . The electron temperature T e 
and emission measure EM are found using the Geometric mean Emission Mea- 
sure (GEM) method. The EM derived from the elements Fe, Si, S, and Al are 
compared in order to measure relative coronal-to-photospheric abundance en- 
hancement factors. We also studied the ion temperature Tj and the non-thermal 
velocity v nt using the line profiles. All these measurements are compared to po- 
lar coronal hole observations from the previous (1996-1997) solar minimum and 
to model predictions for relative abundances. There are many similarities in 
the physical properties of the polar coronal holes between the two minima at 
these low heights. We find that T e , n e , and T\ are comparable in both minima. 
T e shows a comparable gradient with height. Both minima show a decreasing 
T; with increasing charge-to- mass ratio q/M. A previously observed upturn of 
T; for ions above q/M > 0.25 was not found here. We also compared relative 
coronal-to-photospheric elemental abundance enhancement factors for a number 
of elements. These ratios were ~ 1 for both the low first ionization potential 
(FIP) elements Si and Al and the marginally high FIP element S relative to the 
low FIP element Fe, as is expected based on earlier observations and models for 
a polar coronal hole. These results are consistent with no FIP effect in a polar 
coronal hole. 
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Introduction 



A coronal hole is a relatively cool, low-density, and open field-line region of the solar 
atmosphere. During solar minimum, large coronal holes are found at the Sun's polar regions 



as was seen in the 1996-1997 and 2007-2009 solar minima. Coronal 



roles are also known 



source regions of the fast solar wind (IKrieger et al. 



the basic equilibrium form of the solar wind (IBame et al. 



1973 



Zirker 



19771 ). 



19771 ) . which is considered 



There are a number of differences between the recent 2007-2009 solar minimum and 
the previous minimum in 1996-1997. The polar magnetic field during the 2007-2009 solar 



minimum was about 40% weaker and the polar corona 



than during the 1996-1997 minimum (I Wang et al 



hole area was about 20% smaller 



20091 ). Measurements of the solar wind 



from the Ulysses spacecraft show that the fast solar wind during t he 2007-2009 minimu m 



was 3% slower and 1 7% less dense compared to the previous one ( llssautier et al. 



2008 



McComas et al. 



20081 ). Those measurements also show that the solar wind power decreased 



by 25% relative to the 1996-1997 solar minimum. 

In this paper we investigate the physical properties of a polar coronal hole in 2007 at 
low heights and compare them to those of the previous minimum. We ex amine spectral da ta 



from the Extreme U 



2006 



Culhane et al. 



travi olet Imaging Spectrometer (EI S) instrument 



20071 ) onboard the Hinode satellite ( jKosugi et al. 



(IKorendvke et al. 



20071 ). From the EIS 



spectral data we measure the electron density n e , electron temperature T e , emission measure 
EM, and ion temperature T{. We also looked for any "FIP effect". Elements with a low first 
ionization potential (FIP < 10 eV) can undergo what has been dubbed as the FIP effect in 
which the elemental abundances of low FIP ions are enhanced in the corona relative to the 
photosphere. These ratios of coronal-to-photospheric abundances are called "FIP factors" . 
We investigated the FIP effect by comparing ratios of FIP factors. All these results can be 
compared to similar measurements from extreme ultraviolet spectra collected during the 
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1996-1997 solar minimum with the Solar Ultraviolet Measurement of Emitted Radiation 
Spectr ometer (SUMER; 



(CDS; 



Harrison et al 



spacecraft ( 



Wilhelm et al. 



Domingo et al 



19951 ) and the Coronal Diagnostic Spectrometer 



1995 ), bo th on the Solar and Heliospheric Observatory (SOHO) 



19951 ). 



The rest of this paper is organized as follows. The observations are described in 
Section [2j Section [3] reviews the emission measure analysis, which serves as one of our 
primary diagnostics. Data reduction, including line selection and fitting, are described in 
Section HJ Results are presented in Section [5j followed by a discussion and summary in 
Sections El Uncertainties are quoted throughout at an estimated la statistical accuracy. 



2. Observations 



The observations were carried out with the EIS instrument onboard the Hinode 
satellite. The spectrometer covers the wavelength ranges 171 - 211 A and 245 - 291 A and 
has a spectral resolution of 0.022 A per pixel. There are 1", 2", 40", and 266" slit widths 



length of 51 2 ". T 



re observations described here were all 



Brown et al. 



(120081 ) has measured the instrumental line width 



available, each observing a 
performed with the 2" slit. 
(FWHM) for the 1" slit to be 0.054 A in the short wavelength band and 0.057 A in the long 
wavelength band. The instrumental width for the 2" slit is estimated to be 0.007 A greater, 
based on co mparisons betw een observations of the same quiet Sun region with both the 1" 



and 2" slits (jYoung 



2009af ). 



For this study we analyze five observations made on 2007 November 16 at times 06:11, 
06:49, 07:26, 09:40, and 10:50 UT. In each observation the 2"x 512" slit was rastered across 
7 positions in the horizontal direction giving a 14" x 512" field of view. The centers of 
the horizontal scans were at X = —7", 108", 223", 324", and 423" respectively, relative 
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to the center of the Sun (Figured]). Throughout this paper the observations are labelled 
according to their solar- X location. In the vertical direction all the observations were 
centered at Y — 879" relative to the center of the Sun. The line of sight of the portion 
of —7" observation above the limb in the north polar coronal hole is roughly over the 
central meridian. The observations centered at 108" and 223" are further from the central 
meridian, but the above-limb slit portions are within the polar coronal hole. The 324" and 
423" observations cover the quiet Sun corona at low radii and appear to cross into the polar 
coronal hole in the upper portion of the observation. 

For the analysis here we were only interested in the coronal plasma above the limb. 
We therefore chose only pixels above the limb for analysis. As a result the portion of the 
EIS observations used in this study included each of the 7 horizontal pixels (spanning 14" 
at 2" per pixel) and above the limb about 150-250 vertical pixels (1" per pixel). The spatial 
part of the observations below the limb was used to estimate the scattered light intensity 
as described below, but was not considered for the rest of the analysis. 



3. Emission Measure Analysis 



Bryans et al 



fl2009[ ) 



We use the Geometric mean Emission Measure (GEM) method of 
to determine T e and EM based on the intensity of observed spectral lines. The flux of a 
spectral line emitted by a transition from level j to level i of charge state q for element X 
observed from distance d is given by 

1 



And? 



G(T e ,n e )n z c dV. 



(1) 



G(T e ,n e ) is the contribution function defined as 

rijiX+i) n(X +q ) n(X) n(H) A j{ 



G(T e ,n e ) 



n(X +< ?) n(X) ra(H) n e n e 



(2) 
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Here nj(X +q )/n(X +q ) is the relative population of the upper level j for ion X +9 , n(X +q )/n(X.) 
is the relative abundance of charge state q for element X, n(X)/n(H) is the abundance of X 
relative to hydrogen, n(H)/n e is the hydrogen abundance relative to free electrons, and Aji 
is the Einstein A-rate. 

For an isothermal plasma of uniform density Eq. (pQ) simplifies to 

/* = J^G(T e ,rie)EM, (3) 

where EM is the emission measure defined as 

EM = / n 2 c dV. (4) 



The above equations can be solved for the EM as a function of T e and n e giving 

EM = 4n£—^- (5) 

The observed data provide the intensities of the various spectral lines. The function 
G(T c ,n c ) is calculated from atomic data. For most lines G(T e ,n e ) is only a very weak 
function of n e . The EM as a function of temperature can then be plotted for each emission 
line using equation ([5]) where n e is fixed. The density used in the analysis was estimated 
using line intensity ratios (described in Section 15. II) . According to the definition of 
equation (Tjf]) the emission measure depends only on the electron density and the volume 
of the emitting plasma. The assumptions made here imply that these properties should 
be the same for all the ions in the volume, therefore an analysis using the intensity of 
different emission lines originating from the same volume should yield the same EM. The 
temperature is determined by finding the point (EM, T e ) which best satisfies the condition 
that the EM should be the same for each line. When the EM(T e ) curves are plotted, this 
condition is met at the point where the curves intersect. 

The GEM method is a particular way of extracting the averag e values of EM and T e 



using the intersection of the EM(T e ) curves for the observed lines (IBryans et al. 



20091 ) 
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Because EM(T C ) varies rapidly, the GEM approach uses the mean of the logarithm (the 
geometric mean) rather than the arithmetic mean. When multiple intersections occur for a 
given pair of EM curves we eliminate the less physically probable crossings by choosing the 
intersection closest to the average crossing point for all the EM curves. We also eliminate 
crossings between EM curves from the same ion as such crossings are most likely due to 
observational issues or uncertainties in the atomic data. Finally, we omit intersection 
points from temperatures where n(X. +q )/n(X.) is below 1% since under that condition the 



collis ional ionization equilibrium calculations are known to be less reliable ( IBryans et al 
2006k 



G(T e ,n e ) is calculated using the chianti atom ic database ( 



Dere et al. 



Bryans et al. 



1997 



20091) 



fl2009h . The elemental 



and the collisional ionization equilibrium data from 
abundances n(X)/n(H) are assumed to be photospheric, but relative coronal abundances 
can be inferred with further analysis as discussed in Section [5j The photospheric abundances 



are taken from 



Asplund et al. 



( 120091 ) 



4. Analysis of Line Intensities 

The data were reduced using the standard EIS preparation routines to remove the dark 
cu rrent, cosmic r ay spikes, and warm pixels. For warm pixels we followed the procedure 
of 



Youngl ( I2009bl ). The output of the EIS preparation routines converts measured counts 
to physical units (erg cm -2 s _1 sr _1 A -1 ) with an uncertainty based on Poisson counting 
statistics, dark current, and removal of the warm pixels. 

Because emission from the polar coronal hole is faint, to improve statistics we collapsed 
these reduced data into bins consisting of 10 vertical pixels and all 7 horizontal pixels. The 
averaging routine used to collapse the data was careful to ignore bad pixels flagged by the 
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preparation routines. 



Lines were cho s en for a nalysis base d on the identification s of 



Brown et al. 



d2008|), lYoungj (12009d l. and lLandi k YouneJ (120091 ) . Table Q] lists all lines used. 



Young et al. 



(120071 ). 



Note that from here on all wavelengths are given in units of A. We chose lines that were 
either free from significant blends, or where blends were separated in wavelength so that 
the component peaks could clearly be distinguished, or where the blend contribution could 
be subtracted out based on a measurement of a third unblended line known to have a fixed 
intensity ratio with one of the components of the blended line. We also used the Fe x 
A257.26, Fe XII A195.12, and Fe xm A203.82 lines, which are inseparable self blends. In 
each of these cases, since the blend comes from the same ion the total intensity can be used 
in the GEM analysis by summing the contribution functions G(T e ,n e ) for each component 
without introducing additional uncertainties from differences in elemental or charge state 
abundances. 



For each observation we used 25 lines from Fe vm-xill and 8 lines from Si VII, IX, and 
x. In the 423" observation, our furthest outward slit position, the electron temperature 
was greater than at other positions and some Fe xiv lines were strong enough to be used 
in the analysis. For elements other than iron and silicon there were very few lines in these 
observations that were both reasonably strong and unblended. We found 2 S vm lines, 2 
S x lines, 1 Al vm line, 3 Al IX lines, and 3 Mg VII lines. All of the Al lines were very 
weak and could only be analyzed for the lowest few spatial bins above the limb. Both 
S vm lines are blended with Fe XI lines. The Fe XI A198.55 intensity, which blends with 
the S vm A198.55 line, can be estimated by measuring the intensity of the Fe XI A189.13 
line. The theoretical intensity ratio is Fe XI A198.55/A189.13 = 0.84. Similarly, the blended 
contribution from Fe XI A202.63 to the S Vlll A202.61 line can be subtracted by measuring 
the Fe XI A188.23 line intensity where the theoretical ratio Fe XI A202.63/A188.23 = 0.017. 
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Lastly, the Mg VII A278.40 line is a blend with Si VII A278.45, which has a theoretical 
intensity 0.323 times that of the Si VII A275.36 line. In all three cases the pair of lines in 
the ratio originate from the same upper level and hence the intensity ratio is determined 
solely by the branching ratio for transitions to the lower levels. 

The intensities and widths of the emission lines were determined by fitting a sum of 
Gaussian functions plus a linear background to a small wavelength range around each 
line. The Gaussian fits used the eis-auto-fit_gen routine from the EIS analysis software. 
When fitting multiple lines over the same wavelength interval, constraints were included 
as needed in the fitting to keep the position of the line centroid close to the expected 
wavelengths. These constraints typically allowed less than 0.1 A freedom in the centroid 
position, corresponding to a line of sight velocity < 150 kms -1 . This allowed Doppler shifts 
up to this velocity to be detectable. However, we found that centroid positions were very 
close to the expected wavelengths. Comparing fits using slightly different sets of constraints 
confirmed that systematic errors arising from the fitting method did not affect the GEM 
analysis results. Examples of two of the most complex fits are shown in Figure El but the 
majority of fits were of more widely separated lines. 

Line intensities beyond the solar limb consist of coronal line emission plus a component 
from the solar disk scattered within the instrument. Recently the scattered light intensity 
was measured during a partial solar eclipse. It was found that pixels pointed at the eclipsed 
portion of the solar disk read 2% of the int ensity of the pixels observed in the uneclipsed 



portion of the solar disk (jUgarte-Urra 



2010). This scattered light component was removed 
in the current analysis by measuring the intensity of each line for several bins below the 
limb, which was nearly constant, and subtracting 2% of that average intensity from the 
total intensity of the above-limb pixels used in the analysis. For some lines the scattered 
light subtraction brought the intensity in the highest spatial bins to essentially zero. 
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Results 



5.1. Density 

We used the ratio of emission lines from the same ion to measure the electron density. 
Our observations are of a low density region close to the limb and under these circumstances 
photo-excitation and stimulated emissio n from photospheric black body radiation can be 



important processes (j Young et al. 



20031 ). Hence, it was necessary to include these effects in 
the density analysis. Not accounting for these effects gave a density that is too large by 
about 50%. 

In the polar coronal hole observations between —7" to 324" the Fe Vlll A186.60/A185.21 
and Fe IX A189.94/A188.49 intensity ratios show that the density falls from « 8 x 10 7 cm -3 
to ~ 1 x 10 7 cm -3 between 1.05 R & and 1.15 R & (Figure [3]). The density inferred from 
these lines is about the same as what was foun d in polar coronal ho l e measurements from 



the previous 



Fludra et al. 



solar minimum at similar heights ( jBanerjee et al. 



1999 



Landi 



2008|). 



1998 



Wilhelm et al. 



1998 



The Si x A258.38/A261.04, Fe xm (A203.82+A203.79)/A202.04, and Mg VII A280.75/276.15 
ratios, however, yield a density that is roughly a factor of 2 larger at 1.05 R@ and fall off 
slower with height compared to the Fe Vlll and Fe IX density diagnostics. This apparent 
discrepancy for the Si x and Fe xm diagnostics is probably due to the presence of quiet 
Sun coronal plasma along the line-of-sight. For example, assuming the Fe vm and IX lines 
give the actual coronal hole density and assuming typical quiet Sun coronal temperature 
and densities of logT e = 6.15 and n e = 2 x 10 8 cm -3 at 1.05 R@, we estimate that the 
quiet Sun emission measure needs to be only a few percent of the coronal hole emission 
measure in order to produce the inferred densities. This is because the high formation 
temperatures of the Si x and Fe xm lines causes the intervening quiet Sun corona to 
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contribute significantly to the observed intensity. The Fe Vlll and Fe IX ions have a lower 
temperature of formation and are largely insensitive to quiet Sun corona, making them more 
reliable density diagnostics in the relatively cool polar coronal hole. An observational effect 
like this cannot explain the high density given by the Mg VII lines since that ion is abundant 
at typical coronal hole temperatures. We note that discrepancy between the density de rived 



from the Fe vm and Mg VII lines has been observed previously (jYoung fe Landill2009l ). 



Figure [T] shows that the observation at 423" appears to observe the quiet Sun corona 
at low heights. The observation also has a significantly higher temperature, as will be 
discussed in Section 15.21 Consequently, the Si x and Fe xill intensity ratios are the best 
density diagnostics to use in this case. The density in the 423" observation, based on these 
lines, was 2 x 10 8 cm -3 at 1.05 R@ and dropped to about 5 x 10 7 c m" 3 at 1.15 R^. This is 



in agreement with quiet Sun density measurements, such as those of 



Feldman et al 



f|l999h . 



5.2. Electron Temperature 

To measure the temperature as a function of height we performed the GEM analysis 
using the selected iron and silicon for each observation (see Tabled]). These elements had 
several charge states and a large number of lines available so we consider them to be more 
reliable indicators of the temperature than sulfur, aluminum, or magnesium for which there 
were few observed lines. The GEM analysis was performed for each of the 7 pixel x 10 pixel 
spatial bins that were above the limb brightening, usually at distances from the center of 
the Sun R > 1.03 R Q . The density parameter for G(T e ,n e ) versus height was set using n e 
from Fe Vlll and IX for the —7" to 324" observations and using Si x and Fe xill for the 
423" observation. We also performed the analysis with input densities that were larger and 
smaller by an order of magnitude and found that while a small number of lines are density 
sensitive, the average EM and T e are insensitive to the input density. 
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Examples of the EM analysis are shown in Figures H] and [5] for iron and silicon, 
respectively. The EM(T e ) curves for lines of each charge state are in reasonable agreement. 
Measuring temperature in units of K here and throughout the paper, the crossing points 
also appear in a tight cluster near logT e = 6, which supports the idea that the plasma is 
approximately isothermal. The iron line analysis shows that there are some crossing points 
that extend to higher temperatures, which could indicate a higher temperature component 
or a systematic uncertainty. 

Figures E] and [7] show the temperature profiles from iron and silicon, respectively. The 
temperature measurements using either the Fe or the Si lines were in agreement within the 
uncertainties, but the iron T e is systematically greater than that from silicon (e.g., Figure [8]). 
The iron lines come from a larger number of charge states having temperatures of maximum 
abundance in the range log T e = 5.6 — 6.3. In contrast, the silicon lines come from only three 
charge states with temperatures of maximum abundance in the range logT e = 5.8 — 6.1. 
The result is that there is a larger spread in the EM(T C ) crossing points for Fe than for Si 
and a correspondingly greater uncertainty. The large number of iron lines with a maximum 
abundance logT e > 6.1 appear to systematically increase the average temperature derived 
from the iron lines compared to the silicon lines. Fe xn-xiv have maximum abundances at 
log T > 6.2. Removing these lines from the analysis removes the systematic temperature 
difference. 

The three observations in the polar coronal hole at —7", 108", and 223" all showed 
similar temperature profiles, increasing from about logT c ps* 5.95 near 1.03 R Q to 
logT c ps 6.02 near 1.15 R Q (Figures |6] and [7J). The temperature increases steadily with a 
gradient A logT e /Ai? ps 0.75 Rq 1 up to about 1.12 R Q and then levels off above that height. 



These inferred temperatures are very simi 



minimum (IDoschek et al 



1998 



Feldman et al- 



ar to measurements from the previous so lar 



1998 



Fludra et al. 



1999 



Doschek et al 



13 



2 Q ll ; lLandill2008l ). Earlier polar coronal hole observations also showed that the t empera ture 



i ncrea ses with height. We find a temperature gradient very similar to that of 



Landi 



( 120081 ). A possible systematic cause of the apparent temperature increase is the presence of 
surrounding hot quiet Sun coronal plasma along the line-of-sight. Since the scale height 
of the higher temperature plasma is larger, the intensity from intervening hot, quiet Sun 



plasma would decrease with height more slowly tha n the colder polar coron a. 



his w ould be observed as a temperature gradient ( IFeldman fc Landi 



20081 ). 



hole plasma. 



Doschek et al. 



( 120011 ) performed a DEM analysis for a similar observation taken during the previous solar 
minimum and also found a temperature gradient. Their DEM showed a peak in the EM 
distribution at logT e ~ 5.95, but the ratio of high temperature to low temperature DEM 
increased with height. 

In the 324" observation, T e was roughly constant over the height range of the 
observation with logT e w 6.03 — 6.05. Below about 1.10 R Q , T e is clearly greater than in the 
—7" to 223" polar coronal hole observations, but above 1.10 R & the temperature matches 
that of the —7" to 223" observations. Since 324" is near the edge of the polar coronal hole 
this behavior could be caused by quiet Sun corona being observed simultaneously with the 
polar coronal hole at low heights, leading to the calculation of an intermediate temperature 
in our analysis. 

The 423" observation appears to cover the boundary between the quiet Sun corona and 
the polar coronal hole. At low heights the observation falls in the quiet Sun corona and 
the temperature increases from logT e w 6.08 at 1.03 R Q to logT c w 6.15 at 1.10 R & . The 
temperature is then constant up to ~ 1.15 R & where it begins to decline. Figure [1] shows 
that 1.15 R & is close to the intensity transition between the quiet Sun corona and the polar 
coronal hole. Above 1.18 R® the temperature has dropped to about the level of the other 
polar coronal hole observations. The quiet Sun temperature agrees very well with previous 



observations, such as those reported by 



Feldman et al 



(1999) who found logT e « 6.11. 



5.3. Scale Height Electron Temperature 



In every observation the EM as a function of height decreases exponentially (Figure |9]). 
In each case the emission measures and slopes of EM versus height from Fe and Si were 
in agreement to within the uncertainties. The slope of the EM curves is related to the 
scale height and can be used to estimate the temperature. For an isothermal plasma in 
hydrostatic equilibrium the density is proportional to 



where R is radial distance, g is the surface gravity of the Sun, k-Q is the Boltzmann constant, 
and m is the mean particle mass, about 0.61 mn with mn the hydrogen mass. The emission 
measure is a volume integral of the square of the density. If the spatial bins are small, then 
the density does not vary greatly across the field of view and the volume integral reduces 
to an integral along the line-of-sight. Since the radial distance across the observation is 
small, here about 0.15 R & , the length through the plasma along the line-of-sight will be 
approximately constant over this height. The EM is then related to the radius by 



Using this, the electron temperature can be estimated from the slope of a linear fit to 
ln(EM) versus radial distance (Figure [TU]) . The temperatures for each observation using 
this method are given in Table [2j Performing the analysis over shorter sub-intervals did not 
reveal any systematic change in the slope of EM(i2), so we report the results for a linear 
fit over the full above-limb distance. In the polar coronal hole observations —7", 108", 
and 223" the iron and silicon scale height temperatures are greater than the temperatures 
measured by the emission measure analysis (e.g., Figures [61 and [8]). However, the present 



n e oc exp(— mgR/kBT e ), 



(6) 



EM oc exp(-2mgR/k B T e ). 



(7) 
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scale height temperature measurement is closer to T e measured using EM analysis methods 
than has been found in the past. A similar scale height analysis of a polar coronal hole 
during the 1996-1997 solar minimum found a larger difference between the two methods 
with log T e 6.5 from the scale height compared to logT e m 6.0 from an emission measure 



analysis (lLandi 



20081 ). In the 324" observation T e derived from the scale height is about the 
same as that derived from the GEM analysis and in 423" observation the scale height T e is 
smaller than T e from the GEM analysis. 



5.4. FIP factors 



Elemental abundances in the corona may differ from their photospheric values. It is 
commonly observed that elements with a first ionization potential (FIP) below the hydrogen 
Lyman a energy of 10.199 eV are enhanced in the corona relative to the photosphere. This 
is referred to as the FIP effect and the amount of the enhancement is called the FIP factor. 
Based on previous observations and theoretical models, coronal holes are expected to have 



a small FIP effect wit 



( iFeldman et al. 



1998 



i FIP factors of ~ 1 — 2 compared to ~ 3 — 4 in the quie t Sun corona 



Zurbuchen et al. 



2002 



Bryans et al 



2009; 



Laming 



20091 ) 



The EM for two elements can be used to determine ratios of FIP factors. We have 
calculated the EM using photospheric values for the relative abundance n(X)/n(H). Coronal 
abundances can be enhanced by the FIP factor, which is defined as 

f _ Hx)MH)] corona 

/X " [n(X)/n(H)] photosphere - 1 ) 

One can convert G(T e ,n e ) based on photospheric abundances into a G(T e ,n e ) based on 
coronal abundances by multiplying the photospheric abundance based G(T e , n e ) by fx- 
Since EM oc l/G(T e ,n e ) this means EM C = EM p //x, where the subscripts denote quantities 
inferred using coronal (c) versus photospheric (p) abundances. In the logarithm this relation 
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is 

log(EM c ) = log(EM p ) - log(/ x ) . (9) 

If all the emission lines come from the same volume in the corona and the correct coronal 
abundances are used, then EM C should be independent of the element whose emission lines 
are used to calculate it. We can therefore determine the ratio of FIP enhancement factors 
for two elements fx and f\ by 

log(/ Y //x) = log[EM p (Y)] - log[EM p (X)]. (10) 



5.4-1- Iron and Silicon 

Both Fe and Si are low FIP elements with FIPs of 7.902 eV and 8.152 eV, respectively. 
Since both are low FIP elements, they are expected to have similar FIP factors with small 
differences due to differing ionization and recombination rates. None of the observations 
showed a significant trend in /si//Fe as a function of height. An example showing the 223" 
observation is presented in Figure [HJ The uncertainty weighted mean values (/si//Fe) for 
each observation are tabulated in Table [31 The mean is calculated in the logarithm where 
the uncertainties are symmetric. For all polar coronal hole observations here, the ratio 
was consistent with 1. In the quiet Sun corona in the 423" observation the ratio is sightly 
greater than 1. 



The present results are in agreement wi t h previous polar coron al hole observations 



which find fsi/ fFe ~ 1 (IFeldman et al. 



1998 



Zurbuchen et al. 



2002 



. Observations of 



the quiet Sun corona find that the fs; / fv r ratio remains c lose to 1 (IFeldman et al. 



Zurbuchen et al 



2002 : 



Bryans et al 



2009|). 



1998 



Laming! (120091 ) presents a model in which the 



FIP effect is caused by the ponderomotive force generated by Alfven waves propagating 
through the chromosphere and being transmitted or reflected within coronal loops. 
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This model predicts /si//Fe ~ 1 in a coronal hole and fsi/ fpe = 1-0 — 1.5 in quiet Sun 
regions. Both values are consistent with our results and previous observations, although 
the uncertainties in the current measurements limit our ability to make a more detailed 
comparison with theory. 



5.4.2. Sulfur 

Figure [T2] shows the GEM analysis for the four sulfur lines we observed. The 
temperatures derived from these lines are in agreement with the temperatures from the 
iron and silicon lines in every observation. The emission measures were used to calculate 
/s//Fe- The mean values for (fs/fjfa) are presented in Table El None of the observations 
showed a trend in the ratio as a function of height. In most observations we found 
(/s//rj) ~ 0.91 — 1.39. The upper part of the 423" observation, above the apparent T e 
transition to cooler plasma, had a slightly smaller value, (fs/ fpe) ~ 0.82, but this is within 
the error bars of the measurements for the other pointings. 

As mentioned earlier, there appears to be a small amount of quiet Sun corona along the 
line-of-sight. Since the S x lines are more abundant at quiet Sun coronal temperatures and 
are therefore more sensitive to such plasmas, any quiet Sun coronal plasma will increase the 
observed intensity of the S x lines relative to the S vm lines. This increases both the EM 
and the T e at which the crossing points occur in the GEM analysis and thereby increases 
the inferred fs/ fve- Thus, the observed fs/ fpe is actually an upper limit for the polar 
coronal hole measurement. 

The first ionization potential of sulfur is 10.360 eV, making sulfur a marginally high 
FIP element. Sulfur is therefore expected to undergo littl e FIP enhance ment relative to low 



FIP elements. In a coronal hole the FIP effect model of 



Laming! (120090 predicts that the 
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FIP effect is weak and neither sulfur nor iron experience a significant enhancement, giving 
is ~ /fc ~ 1 and consequently fs/fFe ~ 1- Our polar coronal hole observations find that 
fs/fFe is close to 1, which suggests no FIP effect was present in agreement with models. 



Some p olar coronal ho 



/s//fc < 1 (IFeldman et al. 



e ob servations during the previous solar minimum showed 



the photospheric e 



19981 ). These results might be explained by the revisions to 
emental abundance data. If we use pho tospheric abundances from 



Feldman &: Laming! (120001 ) rather than 



Asplund et al 



(120091 ) we find f R //W ~ 0.6 - 0.7 in 



the polar coronal hole observation, which is similar to that found by iFeldman et al. 



(1998) 



The FIP effect is expected to be stronger in the quiet S un corona. Th erefore 



it is expected that in the quiet Sun corona fs/fFe < 1- The iLaminel (120091) model 



predicts /r//fr ~ 0-5 — 0.9. Observations of solar equatorial regions (IFeldman et al. 



1998 



Bryans et al 



20091 ). as well measurements of the slow solar wind ( IZurbuchen et al. 



20021 ) were consistent with a FIP effect for sulfur. However, we do not observe a FIP 
effect near the polar coronal hole in the quiet Sun portion of the 423" observation where 

(/s//Fe) 



i i O+0.18 
± - iz -0.15- 



5.4-3. Aluminum 

Figure [14] shows our GEM analysis for aluminum lines. Only four Al lines were 
sufficiently intense that they could be used in our analysis and even then for only a few 
of the lowest spatial bins with R < 1.10 Rq. The lines come from only two charge states 
that give only three crossing points for the analysis. The small number of lines and crossing 
points means that our use of the standard deviation is not an accurate representation of the 
uncertainty. In those locations where the GEM analysis could be applied to the aluminum 
lines the temperature was similar to the results from other lines. 
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Al is a low FIP element with an FIP of 5.986 eV. We calculated the relative abundance 
factors for Al compared to Fe for all observations. The values for (/ai/ fpe) are shown in 
Table [3J For the 423" observation the Al lines could only be measured in the lower (quiet 
Sun) part of the slit. In every observation /ai//fc wa s not significantly different from 1. 



This is consistent with expectations since both Al and Fe are low 



expected to undergo similar FIP effects. For example, the iLamingi (120091 ) model predicts 



P ele ments, which are 



/ai//fc = 1-0 — 1.1 in a coronal hole and 1.2 in the quiet Sun corona. 



5.4-4- Magnesium 

Only Mg VII had strong unblended lines in these observations. Since no other charge 
states of Mg could be used to find crossing points, the GEM analysis method could not be 
applied. As a result we used a different method for estimating the EM. 

If the observed volume is isothermal we can measure the EM from the magnesium 
lines by using the temperature derived from the GEM analysis of a different element whose 
emission comes from the same volume. Here we use the T e derived from the silicon lines 
which have a formation temperature closer to that for the magnesium lines than do the 
iron lines. We can then determine the EM for the magnesium lines by EM = EM(T 6i si), 
where by T c gj we mean the electron temperature derived from the silicon lines, not the 
temperature of the silicon ions themselves. In this case the uncertainty in the EM arises 
not only from the uncertainty in the measured intensity, which affects the magnitude of the 
EM vs. T e curves and is always present, but also from the uncertainty on T ej s;. In order 
to account for this new source of uncertainty we performed a Monte-Carlo uncertainty 
analysis. We generated a Gaussian distribution of temperature values with the mean value 
T e gi and the standard deviation equal to the uncertainty. We then calculated EM(T e ) for 
each line at each temperature in the set. The mean and standard deviation of the resulting 
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values give EM(T ej si) and its uncertainty. A plot of the EM(T e ) curves showing the result 
using this method is shown in figure [113 

The EM could not be accurately measured in the quiet Sun corona portion of the 
423" observation. In that case the temperature is so high that the expected relative ion 



abundance of n(Mg 7+ )/n(Mg) is less than 1% and the collisiona 



calculations are known to be particularly uncertain (IBryans et al 



ionization equilibrium 



2006 



20091 ). 



The fwig/ fFe ratio in the polar coronal hole observations was determined using the 
derived EM values. The mean (fug/ fpe) from each observation are presented in Table El 
The value of fyig/ fee i n the —7" to 223" observations was about 1.4. The ratio was 
somewhat larger in the 324" observation. This is likely due to interference from the nearby 
quiet Sun corona, which seems to be particularly strong in the 324" observation but could 
be a factor in the others as well. Intervening quiet sun corona would increase T e ,Si- However, 
the Mg VII lines are insensitive to the hotter quiet Sun corona. Thus, with T 6j si larger than 
T e in the volume where most Mg VII emission actually occurs we would overstate the EM 
as well as the inferred fu s / fee- 



Magnesium i s a low FIP element with FIP 7.646 eV. Both 



Feldman et al. 



Zurbuchen et al. 



(|l998h and 



(120021 ) found that polar coronal hole and fast solar wind FIP factors for 
magnesium were only slightly larger than iron during the 1996-1997 solar minimum, with 
a value fug/ fee ~ 1-1- The ponderom otive force model for the FIP effect predicts a ratio 



< 1.2 in a coronal hole (ILaming 



20091 ). Thus, the value in the present measurement is larger 



than expected, but seems to agree within the uncertainty of the method. 
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5.5. Ion Temperature 



The ion temperatures were estimated fro m the widths of t he spectral lines. The width 

~3) 



of a line is related to the ion temperature by (jSeely et al 



199^ 



AA 



FWHM 



AAL t + 41n(2) 



A 



2k B Ti 2 



1 1/2 



'Hi 



Here AAfwhm is the full width at half maximum (FWHM) of the line, AAi nst is the 
instrumental width, T\ is the ion temperature, M the ion mass, k-Q is the Boltzmann 
constant, c is the speed of light, and v nt accounts for line broadening from non-thermal fluid 
flows, such as turbulence. 



An effective velocity can be defined as 



^cff 



2k B Tj 
M 



(12) 



To determine v e g we analyzed only the lines in Table [T] that were not affected by blending. 
We also do not consider lines at positions where the scattered light contributes greater than 
25% of the total line intensity. The scattered light comes from outside the field of view and 
can be expected to have different emission line profiles than emission from within the field 
of view. The widths of the lines were determined by the Gaussian fitting routine described 
above. We determined for each line a value of v c q given by 



c I AA FWHM AAj nst 



A 



41n(2) 



(13) 



We used the method of 



Tu et al. 



(119981 ) to calculate the ion temperature. This method 
assumes that all the ions have the same non-thermal velocity. Thus, the maximum value 
for the non-thermal velocity that can be consistent with the measured width of every ion is 
the effective velocity v e g of the narrowest observed line with the assumption that for that 
ion T; = 0. The lower bound on the ion temperature is obtained by setting v nt equal to this 
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maximum value in equation ffT2]) . and solving for 7]. The upper bound can be found by 
assuming that v nt = and again using equation f[T2"j) to solve for 7]. 

We found the maximum value for t> nt in all the observations was in the rang e 



20 — 40 kms . For polar coronal holes during the 1996-1997 minimum 



Tu et al. 



1998) 



Landi fc Cranmer 



found the maximum value for v Qt was in the range 33 — 44 km s 1 and 
( 2OO9I ) found a range of 30 — 50 kms -1 . The present measurements are in line with these 



earlier observations of polar coronal holes. 

In all our observations the lower bound on the ion temperature was greater than the 
electron temperature for all but the narrowest lines where the lower bound was set to zero 
by definition. Given these results, we adopted a more restrictive lower bound by assuming 
that for the minimum ion temperature 7] = T e , rather than zero. We then proceeded as 
before using a corresponding upper limit on v nt . With this assumption the upper bound on 
v nt was reduced to « 15 — 35 kms -1 . The ion temperature ranges are plotted as a function 
of charge-to-mass ratio q/M in Figures [TBI and fT71 where q is given in units of e and M in 
amu. The results are qualitatively the same as under the assumption that the minimum 
7] = 0, but the tighter bounds on 2] make the dependence on q/M clearer. 

In the observed north polar coronal hole the relation between 7] and q/M is qualitatively 
the same over the height range of each observation (Figure [TB]) . The temperature of 
low-charged, high-mass ions is log 7] ~ 6.5 — 7.4, about an order of magnitude greater than 
the electron temperature. The maximum 7] for q/m < 0.2 ions here is somewhat larger 
than for the 1996-1997 observations, which found log 7] < 7.0. Both Figures [TBI and [T71 
sh ow that the ion tempera ture is lower for q/M > 0.2. This behavior is similar to that 



of 



Landi fc Cranmerl (120091) for a polar coronal hole during the 1996-1997 solar minimum. 



Landi &: Cranmer! ( 120091 ) found that 7] increases for q/M > 0.25, but such behavior was 



not seen in our analysis. 
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For a given height, the behavior of 7] was similar across the observations (Figure [TTj) . 
In particular there was no significant difference in the portion of the 423" observation below 
« 1.15 R & , where we found a higher T e compared to positions within the polar coronal hole. 

Figure [TBI shows how the line widths (i.e. v e g) increase with height, for a selection of 
lines from ions with different q/M. The rate of increase is similar for each of the ions and 
an increase is found in all observations. There are several potential interpretations for the 
increasing line width. One possibility is that it reflects the increasing ion temperature with 
height as suggested in figure [16j This may be due to more heating. The upper bound of 
T[ increases with height by about A log 7] w 0.2 over the range of these observations. The 
increasing width could also reflect an increasing non-thermal velocity. An increasing v nt 
could be caused, for example, by undamped Alfven waves whose velocity amplitude must 
increase as the density decreases in order to conserve wave energy flux F which is given by 




f (Sv 2 ) B. 

4:71 



(14) 



Here, B is the magnetic field strength, p « n e m p is the mass density with m p the mass of a 



proton, and (Sv 2 ) = 2v 2 t is the mean square velocity perturbation from the wave (IJacques 



1977 



Moran 



20011 ) . Similar in creasing line widths were also observed in coronal holes du ring 



the p revious solar minimum (IBanerjee et al. 
20011 ). 



1998 



Wilhelm et al. 



1998 



Doschek et al. 



6. Discussion and Summary 



These measurements of a north polar coronal hole at low heights from the recent solar 
minimum show broad similarities with polar coronal hole observations from the 1996-1997 
minimum. The results are summarized in Table HI Similarities i n the plasma p r operti es of 



polar coronal holes between the two minima were also found by 



Miralles et al 



(I2OIOI ) in 
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polar co ronal holes above 1.7 R & using the SOHO Ultraviolet Coronagraph Spectrometer 



(UVCS; 



Kohl et al. 



19951 ) 



The measured electron density is roughly the same as that measured for polar coronal 
holes during the previous minimum. In situ observations of the solar wind above 1.4 AU 
show that the density of the fast solar wind decreased by rj 20% compared to the 1996-1997 
minimum. The present coronal observation that the density is comparable between the two 
minima at low heights is not inconsistent with the solar wind measurements since density 
flux can be conserved through changes in the magnetic field. 

The electron temperatures measured using the crossings of the EM curves were in 
excellent agreement in both magnitude and gradient with height above the solar limb 
between the 1996-1997 solar minimum measurements and the present results. The 423" 
observation showed a transition in T e from the quiet Sun corona with logT « 6.15 to 
the polar coronal hole with logT e « 6.03 near 1.15 R@. However, we did not observe any 
significant transitions for other measured properties at that location. Using the emission 
measure decay to estimate the scale height in each observation, we found the present 
observations are closer to hydrostatic equilibrium than was seen for 1996-1997. 

We also compared FIP factors relative to iron and found results for silicon and 
aluminum similar in the polar coronal hole to what was found in earlier solar minimum 
polar coronal hole observations. All three elements have low FIPs and are expected to have 
FIP factors that behave approximately the same way. For the low FIP element magnesium, 
we found a slightly larger FIP factor relative to iron, but the magnesium analysis is subject 
to some additional systematic errors. 

For the moderately high FIP element sulfur, we found fs/fFe ~ 1 m the polar coronal 
hole. This measurement for elements with significantly different FIPs suggests that a FIP 
effect was not present in the polar coronal hole because in the absence of a FIP effect 



-25 - 



we expect fs/fFe ~ 1.0. Previous obser v ations have not found a significant FIP effect in 



coronal holes. However, 



Feldman et al 



( 119981 ) studied a polar coronal hole during the 



previous solar minimum and reported a FIP factor for sulfur fs/fpc < 1 consistent with 
the FIP effect. But they found that no other elements of either low or high FIP, except for 
sulfur, showed any indication of a FI P effect. This may be explained by revisions in the 



m ost recent elemental abund ances by lAsplund et al. 
Feldman Sz Laming ( 200ol ). 



as 



(120091 ) compared to earlier data such 



It would be desireable to perform these FIP factor ratio measurements with a larger 
variety of elements in the near future. We measure ratios of FIP factors relative to iron, 
because the EIS spectrum had a large number of strong iron lines which provided a robust 
measurement. It would be better to measure FIP factor ratios by comparing to high FIP 
elements such such as oxygen, nitrogen, argon, or neon, which are expected to experience 
only a very weak FIP effect. Unfortunately, lines from these elements were not available. 
The EIS spectrum should include lines from the oxygen ions O IV— O VI, but the formation 
temperature for these ions is less than logT e = 5.5, well below coronal temperatures. 
Similarly, argon ions Ar XI and Ar XIV may be found in the spectrum; but at logT e = 6.3 
and 6.5, respectively, the formation temperatures for these ions are much higher than the 
observed temperature. 

The inferred ion temperature decreases with increasing ion q/M up to the maximum 
observed q/M < 0.32. For q/M < 0.25 the dependence of Tj on q/M is similar to what was 
observed in the previous solar minimum. Measurements from the 1996-1997 solar minimum 
showed an upturn in 7] for q/M > 0.25, in contrast to the present results where no upturn 
was observed. In the polar coronal hole we found a maximum ion temperature of log Tj ~ 7.4 
for ions with charge to mass ratios q/M < 0.2. The upper bound is slightly higher than 
in measurements from the previous solar minimum, but the inferred temperature range is 
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large so that the possibility of 7] being comparable in the 1996-1997 and 2007-2009 minima 
cannot be ruled out. 

In each observation the line widths were found to increase with height. This can be 
attributed to more ion heating, less cooling, or an increase of the non-thermal velocity. 
In the present observation the height range was limited by the actual field of view of the 
detector and by the low intensity and scattered light contribution at high altitudes. It would 
be interesting to observe the behavior of T\ at higher heights, which requires coronagraphic 
occultation such as is done by UVCS. 

It is surprising that the solar minimum observations are so similar given that polar 
coronal holes are known sources of the fast solar wind and significant changes have been 
observed in the solar wind power, density, and velocity compared to the previous minimum. 
Such changes were not observed in our measurements of the lower regions of the north polar 
coronal hole. 

We thank S. R. Cranmer and J. M. Laming for stimulating discussions. MH and DWS 
were supported in part by the NASA Solar Heliospheric Physics program. The work of EL 
and MPM is supported by several NASA grants. The work of PB was performed under 
contract with the Naval Research Laboratory and was funded by NASA. 
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Table 1:: Line List 



Ion 


A (A) 1 


Transition 1 


Mg VII 


276.154 


o„2 9 2 3 p n 
Zb Zp iq 


9 s 9n 3 3 

za zjl> 01 


Mg VII 2 


278.404 


Zb Zp ±2 


Zo Zp D\ 


Met vn 

o 


280.742 


9 2 o 2 1 7-) 
Zb Zp LJ2 


9 s 9n 3 1 P, 

za zjy rj 


Al viii 


250.139 


9n 2 3 P„ 
Zp 


— 9 c 9r> 3 3 <7, 
Zo Zj7 


Al ix 


280.135 


9 S 2 9« 2 Pi /a 

zb zp r \ 2 


_ c 9-n 2 2 Po ;o 
z* z/; J 3/2 


Al ix 


284.025 


9 s 2 o„ 2 p 
Zi> z/; J3/2 


-9s 9t> 2 2 Po ;o 
Zi> Zp 13/2 


Al ix 


286.376 


9 s 2 9„ 2 p 

za z/; J 3/2 


- 9 s 9t> 2 2 P, /o 
Z5 z/; -t 1/2 


Si vn 


272.648 


9 „2 n 4 3 p 


-9s 9r> 5 3 P n 

Zo Zip 1 \ 


Si vn 


275.361 


2s 2 2 4 3p 


- 9 s 9r> 5 3 Po 

Zj7 1 2 


Si vn 


275.676 


9 q 2 9n 5 3 P, 


-9s 9n 5 3 Pi 

Zb Zp A \ 


Si ix 


258.082 


9o2 9 2 1 7-) 
Zo zp J-^2 


— 9s 9n 3 1 D n 
zb zp AJ2 


Si x 


258.371 


Oq 2 On 2 Po /a 

Zb zp -T3/2 


-9s 9n 2 2 Po /n 
z& z/; J 3/2 


Si x 


261.044 


9 s 2 9 2 p 
z;> z/; J3/2 


-9s 9« 2 2 Pi /n 
Z6 zp j 1/2 


Si x 


272.006 


9,292 p 
zs z/; r 1/2 


- 9s9T, 22C ?i /a 

zb zp J\/2 


Si x 


277.278 


?? 2 ?r> 2 Po /a 

Zb zp J3/2 


-9s 9-n 2 2c v, /a 
z;> z^/ iJl/2 


S viii 3 


198.554 


2s 2 2p 5 2 P 3/2 


- 2s2p b ^i /2 


S VTTT 4 
O V 111 


9D9 fi1 n 


2s 2 2p 52 P 1/2 


-2s2p 62 S 1/2 


S X 


257.147 


2s 2 2p 34 S 3 /2 


-2s2p 44 P 1/2 


S x 


264.231 


2s 2 2^^3/2 


-2 S 2p 44 P 5/2 


Fe viii 


185.213 


3/ 3d 2 L> 5 /2 


-3p 5 3rf 2 ( 3 F) 2 F 7/2 


Fe viii 


186.599 


3p 6 Sd 2 /^ 


- 3p 5 3d 2 ( 3 F) 2 F 5/2 


Fe viii 


194.661 


3/ 3^^5/2 


-3j9 6 4 P 2 P 3/ 2 


Fe viii 


195.972 


3p 6 3^3/2 


-3p 6 4p 2 P 1/2 
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Ion 


A (A) 1 


Transition 


i 


Fe vill 


197.362 


3p 6 3d 2 D 5/2 


- 3p 5 3d 2 


Cs) 2 p 3/2 


Fe IX 


188.497 


3s 2 3p 5 3d 3 F 4 


- 3s 2 3p 4 


{ 3 P) 3d 23 G 5 


Fe IX 


189.941 


3s 2 3p 5 3d 3 F 3 


- 3s 2 3p 4 


{ 3 P) 3d 23 G 4 


Fe IX 


197.862 


3s 2 3p 5 3d 1 P 1 


- 3s 2 3p 5 4p 1 S 


Fe x 


182.307 


3s 2 3p 52 P 1/2 


- 3s 2 3p 4 


( 3 P) 3d 2 P 3/2 


Fe x 


184.537 


3s 2 3p 52 P 3/2 


-3s 2 3p 4 


{ 1 D)3d 2 S 1/2 


Fe x 


190.037 


3s 2 3p 52 P 1/2 


- 3s 2 3p 4 


{ 1 D)3d 2 S 1/2 


Fe x 


193.715 


3s 2 3p 52 P 3/2 


- 3s 2 3p 4 


('S) 3d 2 D 5/2 


Fe x 


r 257.259 


3s 2 3p 52 P 3/2 


- 3s 2 3p 4 


( 3 P) 3d 4 D 5/2 


^ 257.263 


3s 2 3p 52 P 3/2 


- 3s 2 3p 4 


( 3 P) 3d 4 D 7/2 


Fe XI 


180.408 


3s 2 3p 4 3 P 2 


- 3s 2 3p 3 


{ 4 S) 3d 3 D 3 


Fe XI 


182.169 


3s 2 3p 4 3 Pi 


- 3s 2 3p 3 


{ 4 S) 3d 3 D 2 


Fe XI 


188.232 


3s 2 3p 4 3 P 2 


-3s 2 3p 3 


{ 2 D)3d 3 P 2 


Fe XI 5 


188.299 


3s 2 3p 4 3 P 2 


- 3s 2 3p 3 


( 2 D)3d 1 P 1 


Fe XI 


189.719 


3s 2 3p 4 3 P 


- 3s 2 3p 3 


{ 2 D)3d 3 P 1 


Fe XII 


192.394 


3s 2 3p 34 S 3/2 


-3s 2 3p 2 


( 3 P) 3d 4 P 1/2 


Fe XII 


193.509 


3s 2 3p 34 S 3/2 


-3s 2 3p 2 


( 3 P) 3d 4 P 3/2 


Fe XII 


j- 195.119 
1 195.179 


3s 2 3p 34 S 3/2 
3s 2 3p 32 D 3/2 


-3s 2 3p 2 
-3s 2 3p 2 


( 3 P) 3d 4 P 5/2 
[ 1 D)3d 2 D 3/2 


Fe xill 


202.044 


3s 2 3p 23 P 


- 3s 2 3p3d 3 P 1 
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Ion 


A (A) 1 


Transition 1 




r 203.797 


3s 2 3p 23 P 2 - 3s 2 3p3d 3 D 2 


Fe xill < 






*• 203.828 


3s 2 3p 23 P 2 - 3s 2 3p3d 3 D 3 



2009|). 



1 Wavelengths and transitions given by chianti ( IDere et al.l 
2 Si VII A278.44 blend subtracted as described in Section HI 
3 Fe XI A198.55 blend subtracted as described in Section [U 
4 Fe XI A202.6 3 blend subtracted as described in Section HI 
iBrown et al.l (120081 ) give the upper state as 3 Si. 

None of the lines where blends were subtracted (annotated 1-4) nor the self-blended 
lines (marked by a curly bracket) were used in the ion temperature analysis. 
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Table 2:: Electron temperatures based on the scale height 
of logEM(fl). 



Observation logT c Fe lines logT Si lines 



-7" 6.24 ±0.14 6.08 ±0.03 

108" 6.23 ±0.12 6.16 ±0.06 

223" 6.17 ±0.08 6.08 ± 0.02 

324" 6.10 ±0.07 6.01 ±0.02 

423" 6.04 ±0.04 6.00 ± 0.01 
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Table 3:: Weighted mean FIP factor ratios. 



Ratio Observation 

-7" 108" 223" 324" 423" 

(CH) (CH) (CH) (CH) R > 1.15 R & (CH) R < 1.15 R Q (QS) 

(f Si /f Fc ) 1.151™ 1.101°;- 1.221™ 1.241™ 1.451^ 1.421™ 

(/ S // Fe ) i.i7ig : g i.04±g:ig 1.171™ i-03ig : « 0.8218$ 1.121^ 

(/Ai// Fe ) i.38i™ i.66i^ 1.51+S i.i9l8:i ••• 1-4318:1 

(/ Mg // F e) i.32i8;i* 1.5718:18 i.4iigf 7 2.4418;™ i.ssi^ 



(CH) - The measurement lies in the polar coronal hole. (QS) - The high T c at this 
location suggests the measurement may be in the quiet Sun corona. 
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Table 4:: Observed polar coronal hole properties compared to some represen- 
tative values for the 1996-1997 solar minimum. 



Quantity 



2007 a 



1996-1997 



Reference 1 " 



n c (1.05 Rq) 

GEM log T c (1.05 R e ) 

AlogT c /AR 

Scale Height log T c 

f Si/ fee 

/s//Fe 
/Al//Fe 
/Mg//Fe 
fnt 

logTi(g/M < 0.2) 
T ; vs. q/M 
Av cB /AR 



8 x 10 7 cm -3 
5.97±0.02 c 
fa 0.75 R~' 
6.05 



'0 
6.38 



-, 1t .+0.2l d 
1 - lo -0.19 

-, 17 +0.22d 
1 - L < -0.19 



1.38 



+0.71d 
-0.47 



,+0.34d 
-0.27 



1.32I u _ 
< 40 kms" 1 
6.5 - 7.4 
decreasing 
70 kms" 1 R" 1 



7 X 10 7 cm -3 
5.95 ± 0.04 



G 



i 0.73 R 
6.54 

Rj 1 

w 0.7 



< 50 kms^ 1 
6.3 - 7.0 

decreasing <?/M < 0.25; increasing q/M > 0.25 



40 - 70 krns" 1 R, 



Wilhclm ct al. (1998) 
Landi (2008) 
Landi (2008) 
Landi (2008) 
Feldman et al. (1998) 
Feldman et al. ( 1998) 



Feldman ct al. (1998) 
Landi fc Cranmer (2009) 
Landi fc Cranmer (2009) 







Landi fc Cra nmer (2009) 
Wilhelm et al. (1998) 



a This paper 

b 1996-1997 minimum 

c Based on GEM analysis of silicon lines. 

d From the —7" observation. 
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Fig. 1. — S olar and Heliospheric Obse rvatory (SO HO) Extreme ultraviolet Imaging Tele- 



scope (EIT, 



Delaboudiniere et al. 



1995 ) observations at 171 A and 195 A, corresponding to 



lines from Fe IX/X (logT c = 5.9) and Fe XII (logT c = 6.2) lines, respectively. Here T e is 
measured in units of K. The 14" x 512" windows showing the EIS observations are outlined. 
Only the portion of the observations above the solar limb was used in the analysis. The 
observations are labelled in the text by their horizontal position relative to the center of the 
Sun. The line-of-sight of the observation centered at —7" is over the central meridian and 
from left to right the observations are centered at —7", 108", 223", 324", and 423". 
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Fig. 2. — Examples of emission line fits using multiple Gaussians. Shown are Fe XI 189.72 A, 
Fe IX 189.94 A, and Fe x 190.04 A lines (left) and Fe xi 188.22 A, Fe XI 188.30 A, and Fe ix 
188.49 A lines (right), which were used in the analysis. Both figures are for a height R 
~ 1.08 Rq in the 223" observation. These figures represent two of the most complex fits. 
In most cases lines were unblended or were more widely spaced than in the examples shown 
here. 
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Fig. 3. — Density profile for the 223" observation using several different density diagnostic 
line ratios. The Fe Vlll A186.60/A185.21 and Fe IX A189.94/A188.49 ratios yield a density 
in the range expected for polar coronal holes. The Si x A258.38/A261.04, Fe xill (A203.83 + 
A203.80)/A202.04, and Mg VII A280. 75/276.15 diagnostics imply a density a factor of ~ 2 
larger. The discrepancy for the Si x and Fe xill diagnostics is likely an observational effect 
caused by the presence of a small amount of hotter quiet Sun corona along the line-of-sight 
that affects the higher temperature Si x and Fe xill lines but not the cooler Fe vm and IX 
lines. However, Mg VII is abundant at typical coronal hole temperatures and the reason for 
the discrepancy between the Mg VII line and the Fe vm and IX lines is unknown. The error 
bars show only the statistical uncertainty from the fit. Overall uncertainties can be inferred 
from the scatter in the density profiles. 
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Fig. 4. — Emission Measure EM(T e ) for the iron lines at a height of 1.06 R & in the 223" 
observation. The dashed lines show the mean values and the dotted lines show the la 
uncertainty ranges of T e and EM. Here logT e = 5.99 ± 0.05 and log EM = 42.81 ± 0.25. 
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Fig. 5. — Same as Figure H] but showing the EM(T e ) from the silicon lines. Here logT e = 
5.97 ± 0.01 and log EM = 42.84 ± 0.05. 
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Fig. 6. — Electron temperature profiles for the various observations from the GEM analysis 
using iron emission lines. The la uncertainties in logT e are typically ±(0.05 — 0.07). 
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Fig. 7. — Same as Figure EJ but using the silicon lines. The la uncertainties in logT e here 
are typically ±(0.02 - 0.04). 



-40- 



6.3 



6.2 - 




H 6.1 



6.0 - 



5.9 



1.00 



<> 
A 



<> 



<> <> 
A A 



<> 
A 



A 



<> 



<> 
A 



<> <J> 
A 



1.05 1.10 1.15 1.20 

Distance from Sun Center (R ) 



1.25 



Fig. 8. — Temperatures from the GEM analysis of iron and silicon lines for the 223" ob- 
servation. The black and red dashed horizontal lines show the T e inferred from the scale 
height analysis of iron and silicon, respectively. The dotted lines illustrate the associated 
uncertainties on the scale height T e . The scale height temperature results are summarized 
in Table 
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Fig. 9. — Emission measure obtained using the GEM method for the iron lines in each 
observation. The la uncertainties on log EM are typically ±(0.2 — 0.3). 
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Fig. 10. — The emission measure and la uncertainties obtained from the analysis of the 
iron lines in the 223" observation. The linear fit was used to measure the scale height and 
estimate T e . In this case logT e = 6.17 ± 0.08. 
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Fig. 11. — Ratio of the FIP factors /si/iFe plotted as a function of height for the 223" 
observation. The weighted mean and la uncertainty are indicated by the solid and dashed 
horizontal lines, respectively. The values for each observation are tabulated in Table |3] 
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Fig. 12. — Same as Figure H] but showing EM(T e ) from the sulfur lines. Here logT e = 
5.99 ± 0.02 and log EM = 42.90 ± 0.08. 



-45 - 



6 



4 - 



2 - 




1.00 



1.05 1.10 1.15 1.20 

Distance from Sun Center (R ) 



1.25 



Fig. 13. — Same as Figure ITT1 but for fs/fFe- 
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Fig. 14. — Same as Figure H] but showing EM(T e ) from the aluminum lines. Here logT e = 
5.96 ± 0.02 and log EM = 42.86 ± 0.01. 
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Fig. 15. — EM(T e ) from the magnesium lines at 1.06 R & for the 223" observation. The 
EM was found using the temperature from the analysis of the silicon lines. Here assuming 
log T c = 5.97 ± 0.01 we found log EM = 43.00 ± 0.12. 
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Fig. 16. — Ion temperature range as a function of charge-to-mass ratio for several heights in 
the 223" polar coronal hole observation. The behavior for the other pointings was the same, 
with T\ versus q/M unchanging over the height range of the observations. The data for each 
ion have been offset slightly in q/M for clarity. 
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Fig. 17. — Ion temperature range as a function of charge-to-mass ratio for each observation 
at a height 1.05 R@. The data for each ion have been offset slightly in q/M for clarity. 
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Fig. 18. — The effective velocity v e g in the 223" observation using equation ( fl~3"l) . In each 
case the effective velocity increases with height at a similar rate. This could indicate either 
heating of the ions or increasing non-thermal velocity The statistical error bars on v c q are 
generally smaller than the size of the data points. 
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